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Abstract. In spite of an appearance of some impressive recent results in understanding 
of the superfluid turbulence in Hell they fail to evaluate many characteristics of vortex 
tangle needed for both applications and fundamental study. Early we reported the 
Gaussian model of the vortex tangle in superfluid turbulent Hell. That model is just 
trial distribution functional in space of vortex loop configurations constructed on the 
basis of well established properties of vortex tangle. It is designed to calculate various 
averages taken over stochastic vortex loop configurations. In this paper we use this 
model to calculate some important characteristics of the vortex tangle. In particular 
we evaluate the average superfluid mass current J induced by vortices and the average 
energy E associated with the chaotic vortex filament. 



1 Introduction 

The presence of vortex tangle appearing in the superfluid turbulent Hell essen- 
tially changes hydrodynamic properties of the latter (see e.g.rtp Accord- 
ing phenomena are studied in frame of so called Phenomenological Theory (PT) 
pioneered by Vinen Q] and greatly modified by Schwarz [||. The PT describes 
superfluid turbulence (ST) in terms of the total length of vortex lines (per unit of 
volume) or the vortex line density (VLD) C(t) and of the structure parameters of 
the VT. Knowledge of these quantities allows to calculate some of hydrodynamic 
characteristics of superfluid turbulent Hell such as a mutual friction, sound at- 
tenuation etc. Meanwhile there exist many other physical quantities connected 
to distribution of the filaments and their interaction related with other physical 
phenomena which can not be expressed in terms of the PT. The relevant phe- 
nomena should be covered by appropriate stochastic theory of chaotic vortex 
filaments. Of course, the most honest way to develop such theory is to study 
stochastic dynamics of vortex filaments on the base of equations of motion with 
some source of chaos. However due to extremely involved dynamics of vortex 
lines this way seems to be almost hopeless. Thus, a necessity of a developing 
an advanced phenomenological approach appeared. We offer one variant of such 
approach. The main idea and the main strategy are the following. Although the 
phenomenological theory of the superfluid turbulence deals with macroscopical 
characteristics of the vortex tangle, it conveys the rich information concerning 
the instantaneous structure of the vortex tangle. Namely we know that the VT 
consists of the closed loops labelled by Sj(£), uniformly distributed in space and 
having the total length C(t) per unit of volume. From acoustical experiments 
it follows that filaments are distributed in anisotropic manner and quantitative 
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characteristics of this anisotropy can be expressed by some structure parameters 
(see 0], 0, ||]). Beside this usual anisotropy there is more subtle anisotropy 
connected with averaged polarization of the vortex loops. Furthermore there are 
some proofs that the averaged curvature of the vortex lines is proportional to the 
inverse interline space and coefficient of this proportionality (which is of order 
of unit) was obtained in numerical simulations made by Schwarz [4]. 

The master idea of our proposal is to construct a trial distribution function 
(TDF) in the space of the vortex loops of the most general form which satisfies 
to all of the properties of the VT introduced above. We assume that this trial 
distribution function will enable us to calculate any physical quantities due to 
the VT. In the paper we describe typical shape of vortex loop obtained from 
evaluating of the correlation functions. We also calculate the average hydrody- 
namic impulse (or Lamb impulse) 3y in the counterflowing superfluid turbulent 
Hell and the average kinetic energy E associated with the chaotic vortex loop. 

2 Constructing of the trial distribution function 

According to general prescriptions the average of any quantity (B({sj(£,j)})} 
depending on vortex loop configurations is given by 

<£({s,fc)})> = £ B({s^)})V({s^)}). (1) 

{s 3 fe)} 

Here "P({sj(^ )}) is a probability of the vortex tangle to have a particular configu- 
ration {sj(£j)}. Index j distinguishes different loops. The meaning of summation 
over all vortex loop configurations 2{ s -(£•)} m formula (|j) will be clear from 
further presentation . We put the usual in tiie statistical physics supposition that 
all configuration corresponding to the same macroscopic state have equal prob- 
abilities. Thus the probability ^({sj^j)}) for vortex tangle to have a particular 
configuration {s^ ^-)} should be proportional to l/N a u owec i , where N a a owe d is 
the number of allowed configurations, of course infinite 

« jr~ • (2) 

» allowed 

Under term "allowed configurations" N a u owec i we mean only the configurations 
that will lead to the correct values for all average quantities known from experi- 
ment and numerical simulations. Formally it can be expressed as a path integral 
in space of three-dimensional (closed) curves supplemented with some constrains 
connected to properties of the VT. 

Naiiowed « II / v { s j(0} x constrains {s 3 (^)}. (3) 

3 J 

The constrains entering this relation are expressed by delta functions expressing 
fixed properties of the VT. For instance constrain S((s'j(^)) 2 — 1) expresses that 
parameter £ is the arc length. However this condition will lead to not tractable 
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theory. We will use a trick known from the theory of polymer chains (see e.g.S) , 
namely we will relax rigorous condition and change delta function by continuous 
(Gaussian) distribution of the link length with the same value of integral. This 
trick leads to the following expression for number of way: 



Nallowed OC 



n/^i^xe-^r^. ( 4) 



In the same manner we are able to introduce and treat other constrains con- 
nected to the known properties of the VT structure. The detailed calculations 
are exposed in paper of one of the author , now we write down final expression 
for probability of configurations 

Naiiowed oc J V{s(n)} exp (-£{s(k)}) . (5) 

Here s(k) is one-dimensional Fourier transform of variable s(£) and Lagrangian 
£{s(n)} is a quadratic form of the components of the vector variable s(/c) 

£{s( K )} =5>2(KM a/S (K)sg(«). (6) 

In practice to calculate various averages it is convenient to work with the char- 
acteristic (generating) functional (CF) which is defined as a following average: 



wdPiin)}) = (exp [ - J2J2 P "( K K( 

j K^O 

Due to that our Lagrangian |^is a quadratic form (in s(/t)) and, consequently, 
the trial distribution function is a Gaussian one, calculation of the CF can be 
made by accomplishing the full square procedure to give a result 



W({P d (K)}) = exp - J2 E P?(*)iVf (k)P?(-«) ] • (7) 

V 3 

Elements of matrix N^ v {k) are specified from calculation of total length, anisotropy 
coefficient, curvature and polarization. The explicit form of them is written down 
in @. 

Thus we reached the put goal and have written the expression for trial CF 
which, we repeat, enables us to calculate any averaged of the vortex filament 
configuration. For instance calculating some of the correlation functions we are 
able to describe a typical shape of the averaged curve. It is sketched out in Fig.|l] 




/ 1 



Fig. 1. A snapshot of the averaged vortex loop obtained from analysis of the statistical 
properties. Position of the vortex line element is described as Sj(£j), where is arc 
length, s'j(£j) — dSj(^j) /d£j is a tangent vector, unit vector along the vortex line; 
s 'j (Ci) = d? s j(Zi)/d£j m the local curvature vector; vector production s'j(^j) x s"(£j) is 
binormal which is responsible for mutual orientation of the tangent vector and vector 
of curvature. Close (A£ <C R , where R is the mean curvature ) parts of the line are 
separated in 3D space by distance A£. The distant part (R <C A£) are separated in 3-D 
space by distance y/2nRA^ (with correction due to the closeness) . The latter property 
reflects a random walk structure of the vortex loops. As a whole the loop is not isotropic 
having a "pancake" form. In addition it has a total polarization s'j(£j) x s" (£j - )d£i) 
forcing the loop to drift along vector V n and to produce nonzero superfluid mass 
current in z-direction 



3 Hydrodynamic impulse of vortex tangle 



As an one more illustration to the developed theory we discuss hydrodynamic 
impulse of the vortex tangle Jv which is defined as 



IV 



, PsK 

v 2 



(8) 



The quantity Jv is closely related to momentum of fluid (see ||). The averaged 
< Bj(£j)xs'j(£j) > is immediately evaluated by use of CF (0) to give the following 
result: 

pnha v 



'v 



(9) 



Note that the coefficient includes no fitting parameters but only characteristics 
known from the Phenomenological Theory (see ||). Relation (^) shows that the 
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vortex tangle induces the superfluid current directed against the external super- 
fluid current. It should be expected since there is some preferable polarization 
of the vortex loops. In the experiments this additional superfluid current should 
display itself as suppression of the superfluid density. This effect is 3D analog to 
the famous Kosterlitz-Thoulless effect except of that distribution of the vortex 
lines is not calculated but is obtained appealing to the experimental data. 

Since superfluid density enters an expression for second sound velocity, it 
seems attractive to detect it using transverse second sound testing. To do it we 
have firstly to evaluate transverse change of the p s and, secondly, to develop 
the theory to match it to nonstationary case. The general theory asserts that 
while applying a harmonic external second sound field suppression of superfluid 
density becomes the function of frequency ui of the following form: 

\ SV sJtransv 1 + ^ T J 

Here transverse (SJy/S'V^ )t ra nsv is half of the one given by rel. (||). The quantity 
tj is the time of relaxation of the superfluid current Jy. which is to be found 
from dynamical consideration. First, we have to derive dJy/dt with help of 
the equation of motion of the vortex line elements and, second, to evaluate 
various averaged appearing in right-hand side, function we obtain the following 
final result for change of the second sound velocity. Performing all of described 
procedures one obtains that the relative change A112/U2 of the second sound 
velocity is given by 

^ = -/(T)% (11) 

Here the function f(T) is composed of the structure parameters of the vortex 
tangle 

Apkl f a 2 (I - I xx ) 2 



m = ^ : . ^ , (12) 



Decreasing of the second sound velocity in the counterflowing Hell has been re- 
ally observed about two decades ago by Vidal with coauthors M . Let us compare 



our result (11) with the Vidal's experiment. Using the data on the structure pa- 
rameters one obtains that e.g. for the temperatures 1AAK the value of function 
f(T) is about 620 s 2 /cm 4 . Taking the frequency ui — 4.3 rad/s, used in ||, and 
V ns — 2 cm/s one obtains that A112/U2 ~ 4x 10~ 4 , which is very close to the 
observed value. 



4 Energy of vortex tangle 

In this section we calculate the averaged energy of the stochastic vortex loop 
distributed according trial distribution function ^ . The general expression for 
the energy associated with linear vortices can be written as (see e.g. ||) 
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In 3D Fourier space the average energy E (|l3| ) can be rewritten as 



Li Li 



ik(s.(Si) -sj(0)) 



(14) 



JJ k 



Comparing ([li]) and (Q) it is possible to express the energy E in terms of the 
characteristic Functional 



Li Lj 

- 2 r d 3 k [ f „ S 2 W 



2 ^J^^f^JJ " ^ OT>f(&)WPj&) 



(15) 

Here set of P n (£4) m CF VF({P„(£^)}) is again determined with help of the 
^-functions 

Pn(£n) -0, n9 4*,j 1 ° J 

The relation (|l6| ) implies that we have to choose in integrand in exponent of CF 
only points lying in interval from to £j on i -curve and from to £j on j-curve. 
While evaluation of self-energy of the same loop , i = j , one has to distinguish 
points £j , and to put them to be e.g. and Further results concern the case 
of the only loop of length L. Omitting tremendous calculations we write down 
the final answer in the following form: 

47r a 47r V V 71 " / a o 



pn 2 L 



1 2/ 3 j2 | / 2 

(0F-l) 1/27r5/2c 2 3 ' ' ^/ 2 (v^F-l) 1/2 



47T 

where the quantities / (of order of unit) are expressed via the structure param- 
eters of the VT as follows (below/3 = \J I x — I z /I x ) 



h (/3) = y/2 (3-/3 2 )(arcsin (/?)//?), (18) 
h (/?) = (V^P 2 + (2 - P 2 ) arcsin (/?)//?) /^2(3-^) : (19) 
h {(3) = (2 (3 - (3 2 )f 2 (Vl^W- (arcsin /?)//?) //J (20) 
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Let us comment expression (|l7|). The first term in the right-hand side of (O) 
is just the energy of unit of length of a straight vortex filament (see e.g. H) 
multiplied by its length. In this form it is frequently used in theory of superfluid 
turbulence (see e.g. ||) and in other applications. But there are additional terms. 
The third and forth terms appeared from long-range interaction, they are smaller 
then logarithmic ones (about ten percents). The third term is of especial interest. 
It appeared due to polarization of the vortex loop and its presence implies that 
there is some elasticity of the vortex tangle in V ns direction. Results of the 
previous section showed that the VT induces some additional superfluid flow. 
Therefore one can expect that combination of longitudinal elasticity combining 
with inertia of additional will lead to appearing of elastic waves, 3D analog of 
the Tkachenko waves. 

The second one is also logarithmically large. Logarithmic behaviour points 
out that this contribution came from denominator |s(£) — s(£')|. But it was the 
first (local) term which collected contributions from neighbor points along the 
line. Therefore the third term appeared from accidental self-crossing of remote 
(along the line) parts of the vortex filament. The fact that this term is propor- 
tional to L and is of of the first (local) contribution is due to that the line is 
fractal object with Haussdorf Hd dimension equal Hd — 2. According to general 
theory of fractal lines it has an infinite number of self-crossing with cardinal 
number 2Hd — 3 i.e. it is equivalent to line. 

5 Conclusion 

We briefly exposed an essence of Gaussian model of the vortex tangle and give 
several examples how it can be used for evaluation of important physical char- 
acteristics such as induced momentum and energy of interaction. These charac- 
teristics has been discussed early (see e.g. |j]]), however their evaluation has not 
been performed because of lack of a proper theory. We think that these illustra- 
tions convince that Gaussian model can serve as effective tool to study chaotic 
vortex filaments. 

The work was partly funded by Russian Foundation for Basic Research, Grant 
N 99-02-16942 
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